Active OB stars: structure, evolution, mass loss and critical limits 

Proceedings IAU Symposium No. 272, 2011 © 2011 International Astronomical Union 

C. Neiner, G. Wade, G. Meynet & G. Peters, eds. DOI: 00.0000/X000000000000000X 

The rapid magnetic rotator HR7355f 

Th. Rivinius 1 , R.H.D. Townsend 2 , S. Stefl 1 , D. Baade 3 

^SO, Chile; 2 UW Madison, USA; 3 ESO, Germany 

Abstract. For early type magnetic stars slow, at most moderate rotational velocities have been 
considered an observational fact. The detection of a multi-kilogauss magnetic field in the B2Vpn 
star with P « 0.52 d and v sini « 300 km/s has brought down this narrative. We have obtained 
more than 100 high-resolution, high-S/N echelle spectra in 2009. These spectra provide the 
most detailed description of the variability of any He- strong star to date. The circumstellar 
environment is dominated by a rotationally locked magnetosphere out to several stellar radii, 
causing hydrogen emission. The photosphere is characterized by surface chemical abundance 
inhomogeneities, with much stronger amplitudes, at least for helium, than slower rotating stars 
like a Ori E. The highly complex rotational line profile modulations of metal lines are probably 
a consequence the equatorial gravity darkening of HR7355, and thus may offer an independent 
measurement of the von Zeipel parameter f3. 
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■ 1. Introduction 

The Helium-strong star HR7355 has recently been found to host a strong magnetic 
field. The short period of P = 0.52 d puts it among the shortest period non-degenerate 



magnetic stars known, while its vsini w 310km s 1 is the highest for any known non- 



degenerate magnetic star. 



HR 7355 is a star for which effects like gravity darkening and oblate deformation cannot 
be ignored anymore. This means traditional analysis methods to derive stellar parameters 
will, at best, give uncertain results, and at worst misleading ones. 

The presence of photospheric abundance pattern is intriguing, as in a hot star rota- 
tionally induced meridional flows should dilute such pattern quickly. While it has been 
suggested that the magnetic field would inhibit this circulation, HR7355 is the first case 
where this can be tested observationally: not only are there abundance variations across 
the stellar surface, but the amplitude of the equivalent width of the Hel lines is much 
■ larger that in the similar, though less rapidly rotating, Bp star aOriE. 

o 
o 

2. Observations 

Apart from archival and literature data this work is based on high-resolution echelle 
spectra obtained in 2009 with UVES at the 8.2m Kueyen telescope on Cerro Paranal. 
The instrument was used in its DIC2 437/760 setting, which gives a blue spectrum from 
about 375 to 498 nm and a red spectrum from about 570 to 950 nm, with a small gap at 
760 nm. The slit- width was 0.8arcsec, giving a resolving power of about R = 50 000 over 
the entire spectrum. The UVES observations were done on three chunks, in April, July, 
and September 2009. 

Exposure times were between 120 seconds in April and 30 seconds July to September, 
where the shorter exposures were repeated four times. Since the period is short we did 

f Based on observations under ESO programs 081.D-2005, 383.D-0095. 
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Figure 1. The B3- model vs. observed data: Upper panel: UV IUE short- wavelength spectropho- 
tometry and Stromgren visual photometry vs. modeled fluxes. Middle panel: IUE spectrum vs. 
modeled UV-line profiles. Lower panel: All UVES (red) and FORS2 (blue) spectra averaged vs. 
modeled visual line profiles. 



not average shorter exposures, but treat them as taken at distinct phases. In total, we 
have 104 blue and 112 red spectra. The typical S/N of a single spectrum is about 275 in 
the blue region and 290 in the red. 



3. Stellar parameters 

In order to determine the stellar parameters, we made use of spectral synthesis, both 
for line strength and profiles, as well as for absolute fluxes. The code we used is the 
third, re-programmed version of Townsend's (1997) BRUCE and KYLIE suite, which in 
the further we will call "B3" . 

In a first step, the profile of Cll4267, the strongest of all metal lines, was used to obtain 
the projected rotational velocity of vsini = 310±5kms _1 , in good agreement with 
Oksala et al. (2010). Together with the rotational period of P = 0.5214404 d the equa- 
torial radius then becomes ii^ ?equ . sm ^ = 3.19R©. Assuming that a rapidly rotating 
star is constrained by five independent parameters, usually v eq , z, T e ff , M*, i?*, this leaves 
only three of them to be determined, namely the mass, the effective temperature, and 
the inclination. 

The meaning of the effective temperature T e ff in a rapidly rotating star is not 
straightforward, so we note that T e ff here is the uniform black-body temperature, that a 
star of the same surface area would need to have the same full solid angle luminosity as 
the actual gravity darkened star. The T e ff of HR7355 is constrained using de-reddened 
flux-calibrated UV spectra (from IUE) and visual photometry. These data are fitted best 
with T eff = 17 000 K (see Fig. [J). 
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Figure 2. Phased observational data: Photometric data from Oksala et al.(2010), and ASAS 
data from Mikulasek et al.(2010) (uppermost), EWs from UVES (middle), with FORS and 
FEROS W\ measurements of Hel4388, 4713 added, as well as FORS (Rivinius et al. (2010)) 
and ESPaDOns magnetic data (lowermost panel). 



This temperature seems incompatible with the spectral classification of B2. However, 
due to chemical peculiarity the helium lines are much stronger than they would be for 
solar abundances, biasing classifications towards the spectral type with the strongest Hel 
lines, which is B2. 

Since P rot and vsini are well constrained, the inclination effectively determines the 
stellar equatorial radius and the geometric projection onto the line of sight, i.e. the area 
of the star seen. Together with the above derived T e ^ and the Hipparcos distance, the 
flux levels of the SED for 17 000 K constrain the inclination, therefore. For T e ^ = 17000 
models with i = 60° are just so compatible with the flux curve for the closest possible 
distance; better fitting inclinations result in less plausible parameters elsewhere. 

Although this still does not constrain the stellar mass in a straightforward observa- 
tional way, all the other parameters are sufficiently known to leave only a narrow range 
of acceptable evolutionary track masses. This is because we can, at least in first order, 
expect that the luminosity of a star of a given mass does not change with its rotational 
velocity. The model has a luminosity of about L* = 1000 Lq, and such a luminosity is 
only compatible with a mass of about 6 ± 0.5 M , not with a higher one. 
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Figure 3. Phased variations of photospheric Helium and metal lines constructed with 64 bins. 
4. Variability 

For the ephemeris we adopt the epoch given by Rivinius et al. (2010), Xb(MJD) = 
54 940.33, which is the mid-date of an occupation of the star by the magnetospheric 
material, seen in the Balmer line spectroscopy, and as such very well defined. The best 
period proved to be the value published by Oksala et al. (2010), P — 0.5214404 d, who 
could rely on a twice as long time-base. Figure [2] shows previously published data together 
with equivalent width measured in the UVES data, phased with the ephemeris given here. 

The Ha equivalent widths are fully dominated by variations of the circumstellar en- 
vironment, there is no indication for a photospherically intrinsic variation of this or any 
other hydrogen line's strength. The equivalent width of helium lines like Hel4713 and 
4388 are clearly of a double wave character. Although the strongest EW points in both 
half-cycles are on a similar level, the weakest EW values differ by a factor of two. 

In terms of line profiles, the most obvious photospheric variation is that of the Hel lines. 
The variations are very strong, in percentage of the line strength much stronger than in 
other He- strong stars like a Ori E. Figure [2] shows a factor of three to four between the 
weak and strong states of Hel lines. Other than the figure suggests, however, a detailed 
look at the spectra shows that the variation is not as plain as consisting of two enhanced 
polar regions and a depleted belt only, in particular not in the metal lines, which hardly 
show EW variations, but clear line profile variability (see Fig. [3]). 
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5. Magnetic field and magnetosphere 

Constraints on both the inclination i and the obliqueness /? can be obtained from 
the observed properties of the magnetic field. Applying standard diagnostics of stellar 
magnetic fields gives r = B m [ n / B max = cos(/3 + cos(/3 — i) = —0.78. This means that 
j3 + i < 140°, for which i = /3 = 70°. Since we have determined z = 60° above, we 
conclude that both i and /3 are between 60 and 80°, in such a way that /3-\-i is about 130 
to 140°. The field strength at the magnetic poles would then be 11 kG. As will be shown 
in the next paragraph things are probably more complicated than in a pure dipole field, 
but the above values should at be at least a good approximation. 

The measured magnetic curve, when fitted with a sine, has an offset of +0.3 kG, mag- 
netic nulls expected to occur at <p = and = 0.54 (dotted lines in lowermost panel of 
Fig. [2]). Yet, the cloud transits are observed to occur at (j) = 0.0 and <p = 0.46 (dotted lines 
in upper panels of Fig. [2]). Similarly, the magnetic poles pointing towards the observer 
would be expected for = 0.28 and 0.78 in case of a sinusoidal variation (dashed lines in 
lowermost panel of Fig. [2}, but spectroscopically, i.e. by maximal He-enhancement, are 
rather seen at (j) = 0.25 and 0.68 (dashed lines in upper panels of Fig. [2]). Both these off- 
sets could be explained with a misaligned magnetic dipole wrt. to the stellar center, which 
would display in a non-sinusoidal variability curve of the magnetic field. Unfortunately, 
there are not enough magnetic measurements to address this point in detail. 

The two occultations occur at phases (j) = 0.0 and <p = 0.46 (Fig. [4]). The passing 
of the lobe is fast, about A</> = 0.13 from — vsmi to -fvsinz. If the magnetospheric 
lobes reached down to the star, this would rather be Acj) = 0.5. As the crossing time 
decreases quadratically with distance, there is no absorbing material directly above the 
star until about 2 R*. A quarter of a cycle later, the emitting material is seen next to the 
star, and since it is in magnetically bound corotation there is a linear relation between 
velocity and distance from the stellar surface. The emission has an inner edge at about 
600kms _1 « 2vsmi, which again points to an empty region inside 2R*. The outer edge 
of the emission is at a velocity of about 4 x vsini in Ha and 3 x vsmi in Pai4, which 
gives outer geometrical limits for the lobe emission of about 4R* and 3R*, respectively. 

The theoretical line profiles from the B3 modeling are good approximations of the pho- 
tospheric profile, as seen in Fig.[U The residuals (Fig. [4j) represent the clean circumstellar 
emission and it is possible to measure Balmer decrements when the emitting material 
is next to the star. While the values for D54 would be in agreement with logarithmic 
particle densities between 11.7 and 12.5 per cm 3 , the densities from ^34 are somewhat 
higher, at 12.2 to 12.8 per cm 3 . In any case, these values are close to the optically thick 
limit, above which the decrements become independent of density. 



6. Conclusions 

Apart from its high v sini, HR7355 is a rather typical member of the class of He- 
strong stars. Due to its relative proximity and brightness it is an ideal target to study the 
interplay of rapid rotation effects, like gravity darkening and meridional circulation, with 
magnetic effects, like chemical peculiarities and inhibition of the meridional circulation 
and other non-rigid motions. Indeed it seems very unlikely that a relatively strong and 
ordered field as seen in HR 7355 could give rise to so many distinct zones of Helium and 
metal abundances. 

The presented data, more than 100 high-quality spectra filling the phase diagram, 
provide an excellent base to apply techniques like Doppler imaging, rigidly rotating mag- 
netosphere models, and Monte-Carlo radiative transfer models. 
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Figure 4. Colour-coded, phased observational variations of the circumstellar magnet osphere. 
Ha and H/3 are the residuals after subtracting the B3 model; the residual level in the line center 
is unity to better than 1 %. The respective original spectra are overplotted for Ha and H/3. 
Emission clearly seen in residuals, in fact up to KS. Paschen lines show emission clearly as well, 
shown here is Pai4; note that these are not residuals and the different abscissa scale. 128 phase 
bins were used for all panels in order to sample the fast variations around the occultation phases 
in sufficient detail. 
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Discussion 

Groh: Could you comment on the expected angular extension of the Ha and Br7 line 
forming regions and whether those could be resolved by interferometric observations? 

Rivinius: Due to the rigidly rotating magnetosphere, R max = ^max/(^sini), so for both 
regions we expect about 3 to 4 stellar radii. For HR 7355 this is well sub-milliarcsecond, 
but for HR5907 (see Grunhut et al., this volume) at least the spectrally dispersed phase 
signature is probably in reach of AMBER/ VLTI and the CHARA combiners. 

Wade: Due to the peculiarities hydrogen may be non-uniformly distributed in the at- 
mosphere. Could the apparent departures from dipolar field topology be due to an effect 
like this? 

Rivinius: I would be surprised, since a) already the timing of the magnetosphere oc- 
cupations, i.e. without looking at the magnetic field at all (except polarity) requires a 
non-sinusoidal (i.e. non-dipole) field signature and b) also with the ESPaDOnS data, 
derived from metal and He- lines, a similar argument can be made. 



